function y = func_H_m2(q, m, CONSTS)

    k0 = CONSTS.k0;
    a = CONSTS.a;
    
%     hankel_m = besselh(m, 2, k0*a*q);
%     hankel_mp1 = besselh(m+1, 2, k0*a*q);
    hankel_m = H_m2_from_arg(k0*a*q, m);
    hankel_mp1 = H_m2_from_arg(k0*a*q, m+1);
    
    y = (hankel_mp1./(k0*a*q.*hankel_m)) -  (m./((k0*a*q).^2));
end
